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175 (Dated: December 9, 2011) 

Abstract 

We present a search for the standard model Higgs boson production in association with a W 
boson in proton-antiproton colhsions {pp — > W^H lubb) at a center of mass energy of 1.96 TeV. 
The search employs data collected with the CDF II detector which correspond to an integrated 
luminosity of approximately 2.7 fb~^. We recorded this data with two kinds of triggers. The first 
kind required high-pr charged leptons and the second required both missing transverse energy 
and jets. The search selects events consistent with a signature of a single lepton (e^//x^), missing 
transverse energy, and two jets. Jets corresponding to bottom quarks are identified with a secondary 
vertex tagging method and a jet probability tagging method. Kinematic information is fed in an 
artificial neural network to improve discrimination between signal and background. The search 
finds that both the observed number of events and the ncTiral network output distributions are 
consistent with the standard model background expectations, and sets 95% confidence level upper 
limits on the production cross section times branching ratio. The limits are expressed as a ratio 
to the standard model production rate. The limits range from 3.6 (4.3 expected) to 61.1 (43.2 
expected) for Higgs masses from 100 to 150 GeV/c^, respectively. 

176 PACS numbers: 13.85.Rm, 14.80.Bn 
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177 I. INTRODUCTION 



17B Standard electroweak theory predicts the existence of a single fundamental scalar particle, 

179 the Higgs boson, which arises as a result of spontaneous electroweak symmetry breaking 

180 The Higgs boson is the only fundamental standard model particle which has not been exper- 

181 imentally observed. Direct searches at LEP2 and the Tevatron have yielded constraints on 

182 the Higgs boson mass. LEP2 data exclude a Higgs boson with rriH < 114.4 GeV/c^ at 95% 

183 confidence level (C.L.). Recently, the Tevatron has excluded at 95% CL. the mass range 

184 154 < rriH < 175GeV/c^ In addition, recent global fits to electroweak data yielded a 

185 one-sided 95% confidence level upper limit of 158 GeV/c^ jsl. If the experimental lower limit 

186 of 114.4 GeV/c^ is included in the fit, then the upper limit raises to 185 GeV/c^. 

187 The Higgs boson branching ratios depend on the particle's mass. If the Higgs boson has 

188 a low mass {rriH < 135GeV/c^), it decays mostly to bb {4]. If the Higgs boson has a high 

189 mass {niH > 135GeV/c^), then it preferentially decays to W^W~ . 

190 Higgs boson production in association with a W boson (WH) is the most sensitive low- 

191 mass search channel at the Tevatron. WH production is more sensitive than ZH production 

192 because it has a larger cross section. It is more sensitive than direct Higgs production 

193 gg H bb because it has a smaller QCD background. 

194 Searches for WH — )■ iubb at ^/s = 1.96 TeV have been recently reported by CDF using 

195 1.9 fb^^ and DO using 440 pb~^ jc]. The CDF analysis looked for WH production in 

196 charged-lepton-triggered events. It improved on prior results by employing a combination of 

197 different jet flavor identification algorithms j?]. Flavor identification algorithms distinguish 

198 between jets that are induced by light partons {u, d, s, g) and jets containing the debris of 

199 heavy quarks (&, c). The analysis also introduced multivariate techniques that use several 

200 kinematic variables to distinguish signal from background. The analysis set upper limits on 

201 the Higgs boson production rate, deflned as the cross section times branching ratio a ■ B for 

202 mass hypotheses ranging from 110 to 150 GeV/c^. The rate was constrained to be less than 

203 1.0 pb at 95% C.L. for itlh = 110 and less than 1.2 pb for 150GeV/c^. This corresponds 

204 to a limit of 7.5 to 102 times the standard model cross section. More recently, CDF has 

205 produced a search with 2.7 fb~^ of data that combines both neural network and matrix 

206 element techniques [8i]. The search we present here is an ingredient in the most recent 

207 combination. 
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20B The new search for WH — )■ iubb reported here builds on the previous CDF result by 

209 adding more data and introducing new analysis techniques for identifying W candidate 

210 events that have been recorded using triggers involving missing transverse energy I^x and 

211 jets. We use 2.7 fb^^ of data in our search, which is an increase of nearly 50% over the 

212 prior search. Our analysis uses both events recorded with a charged-lepton trigger and 

213 events recorded by a trigger that selects missing transverse energy I^t and two jets. The 

214 missing transverse energy vector is the negative of the vector sum of calorimeter tower 

215 energy deposits in the event. It is corrected for the transverse momentum of any muons in 

216 the event. I^t is the magnitude of the missing transverse energy vector. Missing transverse 

217 energy suggests that a neutrino from a W decay was present in an event. We identify W 
21B candidates in I^t + jet events using looser charged-lepton identification requirements that 

219 recover muons that fell into gaps in the muon system. We show that including these events 

220 significantly increases the search sample and that these new events have a purity that is 

221 comparable to the samples using charged-lepton triggers samples. 

222 We describe the analysis as follows: in Section [III we describe the CDF II detector. We 

223 explain the event selection criteria in Sec. IIIIl focusing especially on the identification of loose 

224 muons. In Sec. IIII Dl we discuss the 6-tagging algorithms. We estimate contributions from the 

225 standard model (SM) backgrounds and show the results in Sec. [IVl In Sec. |Vl we estimate our 

226 signal acceptance and systematic uncertainties. Sec. |VT] describes the multivariate technique 

227 that we use to enhance our discrimination of signal from backgrounds. We report our 

228 measured limits in Sec. IVIII and interpret the result in Sec. IVIIII 

229 II. CDF II DETECTOR 

230 The CDF II detector [9] geometry is described using a cylindrical coordinate system. The 

231 z-axis follows the proton direction, the azimuthal angle is 0, and the polar angle 6 is usually 

232 expressed through the pseudorapidity r] = — ln(tan(^/2)). The detector is approximately 

233 symmetric in rj and about the z axis. The transverse energy is defined as = E sin 9 and 

234 transverse momentum as Pt = psin^. 

235 Charged particles are tracked by a system of silicon microstrip detectors and a large open 

236 cell drift chamber in the region |?7| < 2.0 and |?7| < 1.0, respectively. The open cell drift 

237 chamber is called the central outer tracker (COT). The tracking detectors are immersed 
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238 in a 1.4 T solenoidal magnetic field aligned coaxially with the incoming beams, allowing 

239 measurement of charged particle momentum. 

240 The transverse momentum resolution is measured to be 6pt/pt ~ 0.1% ■ pT(GeV) for 

241 the combined tracking system. The track impact parameter is the distance from the 

242 event vertex to the track's closest approach in the transverse plane. It has a resolution of 

243 o"((io) ~ 40 /im of which 30 /xm is due to the size of the beam spot. 

244 Outside of the tracking systems and the solenoid, segmented calorimeters with projective 

245 tower geometry are used to reconstruct electromagnetic and hadronic showers 10Nl2| over 

246 the pseudorapidity range \ri\ < 3.6. A transverse energy is measured in each calorimeter 

247 tower where is calculated using the measured z position of the event vertex and the tower 

248 location. 

249 Small contiguous groups of calorimeter towers with energy deposits are identified and 

250 summed together into an energy cluster. Jets are identified by summing energies deposited 

251 in electromagnetic (EM) and hadronic calorimeter (HAD) towers that fall within a cone 



252 of radius AR = ^ {A(f)'^ + At]'^) < 0.4 units around a high-ii^T' seed cluster 13|]. Jet en- 

253 ergies are corrected for calorimeter non-linearity, losses in the gaps between towers and 

254 multiple primary interactions [l^. Electron candidates are identified in the central electro- 

255 magnetic calorimeter (GEM) as isolated, electromagnetic clusters that match a track in the 



256 pseudorapidity range \ri\ < 1.1. The electron transverse energy is reconstructed 



257 electromagnetic cluster with a precision a{ET)/ET = 13.5%/ J Er/iGeV) © 2% [lOj 



rom the 



258 This analysis uses three separate muon detectors and the gaps in between the detectors to 

259 identify muon candidates. After at least five hadronic interaction lengths in the calorimeter, 

260 the muons encounter the first set of four layers of planar drift chambers (CMU). After passing 

261 through another 60 cm of steel, the muons reach an additional four layers of planar drift 

262 chambers (CMP). Muons require pt > 1.4GeV/c to reach the CMU 15| and an px > 2.0 

263 GeV/c to reach the CMP jl^. Muon candidates are then identified as tracks that extrapolate 

264 to line segments or "stubs" in one of the muon detectors. A track that is linked to both 

265 CMU and CMP stubs is called a CMUP muon. These two systems cover the same central 

266 pseudorapidity region with |?7| < 0.6. Muons that exit the calorimeters at 0.6 < < 1.0 

267 are detected by the CMX system of four drift layers and are called CMX muons. Tracks 

268 that point to a gap in the CMX or CMUP muon system are called isolated track muon 

269 candidates. 
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270 The CDF trigger system is a three-level filter, with tracking information available even 



271 at the first level 17|]. Events used in this analysis have passed either the electron trigger, 

272 the muon trigger, or the missing transverse energy Efx trigger selection. The lepton trig- 

273 ger selection is identical to the selection used in j5|. The first stage of the central electron 

274 trigger requires a track with > 8 GeV/c pointing to a tower with Ej- > 8 GeV and 

275 E}iAj)/E-EM < 0.125, where -Erad is the hadronic calorimeter energy and Eem is the elec- 

276 tromagnetic calorimeter energy. The first stage of the muon trigger requires a track with 

277 pt > 4: GeV/c (CMUP) or 8 GeV/c (CMX) pointing to a muon stub. For lepton trig- 
27B gers, a complete lepton reconstruction is performed online in the final trigger stage, where 

279 we require Et > 18 GeV/c^ for central electrons (GEM), and pt > 18GeV/c for muons 

280 (GMUP,GMX). 

281 The I^T plus two jets trigger has been previously used in the V{= W,Z)H — )• I^t +bb 

282 Higgs search [l8| and offers a chance to reconstruct WH events that did not fire the high-pr 

283 lepton trigger. The trigger's requirements are two jets and missing transverse energy. The 

284 two jets must have Et > 10 GeV, and one must be in the central region |?7| < 0.9. The 

285 missing transverse energy calculation that is used in the trigger, assumes that primary 

286 vertex of the event is at the center of the detector and does not correct for muons. The 

287 trigger requires I^^^ > 35 GeV. Sections Hn] and |V] discuss the implications of these trigger 

288 requirements on the event selection and trigger efficiency. 



289 III. EVENT SELECTION 



290 The observable final state from WH production and decay consists of a high-p^ lepton, 

291 missing transverse energy, and two jets. This section provides an overview of how we re- 

292 construct and identify each part of the WH decay, focusing especially on isolated track 

293 reconstruction, which is new for this result. Additional details on the event reconstruction 

294 can be found in Ref. Isl. 



295 A. Lepton Identification 

296 We use several different lepton identification algorithms in order to include events from 

297 multiple trigger paths. Each algorithm requires a single high-pr (> 20 GeV/c), isolated 
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298 charged lepton consistent with leptonic W boson decay. We employ the same electron 



igj and 



299 and muon identification algorithms as the CDF W cross section measurement 

300 the prior CDF WH search p\. We classify the leptons according to the sub-detector that 

301 recorded them: CEM electrons, CMUP muons, and CMX muons. We supplement the lepton 

302 identification with an additional category called "isolated tracks". An isolated track event 

303 is required to have a single, energetic track that is isolated from other track activity in 

304 the event and that has not been reconstructed as an electron or a muon using the other 

305 algorithms mentioned above. 

306 The isolated track selection is designed to complement the trigger muon selection in that 

307 it finds muons that did not leave hits in the muon chambers, and therefore, could not have 

308 fired the muon trigger. Figure [T] shows how isolated track events increase overall muon 

309 coverage. The isolated track events are concentrated in the regions where there is no other 

310 muon coverage. Including isolated track events increases the acceptance by 25% relative to 

311 the acceptance of charged-lepton triggers. 



CDF Muon Coverage 




0- 



-3 



CDF Muon Coverage 




-1 .5 



-0.5 



0.5 



1.5 



FIG. 1: (Left) Angular distribution of WH Monte Carlo muon triggered events. Note the cracks 
between CMUP chambers and the gap between the CMUP and CMX. (Right) Isolated track events 
recover high-p^^ muons that fall in the muon chamber gaps. 

We identify isolated tracks based on criteria used in the top lepton plus track cross section 



313 measurement 20|] . Table |T] outlines the specific isolated track selection criteria. The track 



314 isolation variable quantifies the amount of track activity near the lepton candidate. It is 
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Variable 



Cut 



Pt > 20 GeV/c 

l^ol < 60 cm 

Idol < 0.02 cm 

\do\ (no Si hits) < 0.2 cm 

track isolation > 0.9 

Axial COT hits > 24 

Stereo COT Hits > 20 

Num Si Hits (only if num expected hits > 3) > 3 



TABLE I: Isolated track identification requirements. In the table, do is the track impact parameter, 
do (no Si Hits) is the impact parameter for tracks that have no silicon tracker hits, zq is position 
along the direction of the beamline of the closest approach of the track to the beamline, and the 
Axial and Stereo hits are on tracks the open cell drift chamber (COT). We define track isolation 
according to equation [TJ 

315 defined as 

pricandidate) 

316 irklsol = — -. -— — -. — — , (1) 

PTycandtdate) + J2PT{trk) 

317 where J^Pritrk) is the sum of the pt of tracks that meet the requirements in Table HTl Using 

318 this definition, a track with no surrounding activity has an isolation of 1.0. We require track 

319 isolation to be > 0.9. 

320 We veto events with an identified charged lepton that fires the trigger (CEM, CMUP, 

321 CMX) in order to ensure that the data sets are disjoint. In addition, we veto events with 

322 two or more isolated tracks or a single isolated track that falls inside the cone of a jet {AR < 

323 0.4), as these events are unlikely to have come from W fiu decay. 



324 B. Jet Selection 

325 WH signal events have two high-i?T jets from the H ^ bb decays. We define recon- 

326 structed jets using a cone of AR < 0.4, where AR = \/A0M-'A?p. We require jets to have 

327 Et > 20 GeV and |?7| < 2.0. The r] cut ensures that the jets are within the fiducial volume 
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Variable 



Cut 



Pt > 0.5 GeV/c 

Ai2(trk, candidate) < 0.4 

AZ(trk, candidate) < 5 cm 

Number of COT axial hits > 20 
Number of COT stereo hits > 10 



TABLE II: Requirements for tracks included in track isolation calculation. 

Trigger Sample Jet Selection 

Et > 20 GeV 

Charged Leptons 

\r]\ < 2.0 
Et > 25 GeV 
\r]\ < 2.0 

I^T + Jets 

At least one jet \r]\ < 0.9 
AR > 1.0 



TABLE III: Jet selection criteria for events in our different trigger samples. 

328 of the silicon detector. The jet energies are corrected to account for variations in calorimeter 

329 response in t], calorimeter non-linearity, and energy from additional interactions in the same 

330 bunch crossing. Monte Carlo simulations (MC) show that about 60% of WH events passing 

331 our selections result in two-jet events. The remainder is split evenly between events with 

332 one or three jets. Events with one or three jets have a worse signal-to-background ratio than 

333 those with two jets due to contamination from background processes such as VT-l-jets and 

334 ti, respectively. We limit our search for WH — )■ iubb to events with W + exactly two jets. 

335 For events collected on the I^t + jets trigger, we require the jets to have an Et > 25 GeV 

336 to ensure that they are above the trigger threshold. One of the two jets must be in the 

337 central region |?7| < 0.9 to match the requirements of the trigger. In addition, because 
33B the trigger has a low efficiency for jets that are close together, we require the jets to be 

339 well-separated {AR > 1.0). 

340 Table UTTl summarizes the jet selection criteria for events in each trigger sample. 
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341 In calculating event kinematics we find it useful to consider loose jets that have either 

342 somewhat smaller Et than our cuts or have high-ii^T but are further forward than our 

343 standard jets. We call these jets "loose jets" . We do not use them directly in our event 

344 selection, but we do use them in calculating kinematic variables. We define loose jets to 

345 be jets with Ef > 12 GeV in the region |?7| < 2.0, and Et > 20 GeV in the region 

346 2.0 < \r]\ < 2.4. 



347 C. Missing Transverse Energy 

348 The presence of a neutrino from the W decay is inferred from the presence of a significant 

349 amount of missing transverse energy. The missing transverse energy vector is the negative 

350 of the vector sum of all calorimeter tower energy deposits with |?7| < 3.6. The Efrjs, the 

351 magnitude of the missing Et vector. We correct the energy of jets in the event ij] and 

352 propagate the corrections to the EfT ■ We also account for the momentum of any high pt 

353 muons. When we calculate J^t , we use 2;-position of the primary vertex to get the correct Et 

354 for each calorimeter tower. Some events have more than one vertex. In this case. We use the 

355 sum of the transverse momentum of the tracks associated with each vertex to distinguish 

356 between the vertexes. The primary vertex is the one with the highest sum of the track 

357 transverse momentum. We then require I^t to exceed 20 GeV. 



35B D. 6-jet identification 

359 Both of the jets in WH events originate from H ^ bb decays. Many backgrounds have 

360 jets that come from light-flavor partons {u, d, c, s, g), such a.s W + jets and QCD. Jets from 

361 b quarks can be distinguished from light-flavor jets by looking for the decay of long-lived 

362 B hadrons. We use the same 6-jet identification strategy as the previous WH search jsl. 

363 We employ two separate algorithms to identify B hadrons. The secondary vertex tagging 



364 algorithm 2l| takes tracks within a jet and attempts to reconstruct a secondary vertex. If a 



365 vertex is found and it is significantly displaced from the primary vertex, the jet is identified, 



366 or tagged, as a 6 jet. The Jet Probability algorithm j22|] also uses tracking information inside 

367 of jets to identify B decays. Instead of requiring a secondary vertex, the algorithm looks 
36B at the distribution of impact parameters for tracks inside a jet. If the jet has a significant 
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369 number of large impact parameter tracks, then it is tagged as a 6-jet. Jet probability tags 

370 have a lower purity than secondary vertex tags. 



371 E. Lepton + Jets Selection 

372 After identifying the final state objects in the event, we purify the sample with quality 

373 cuts. We fit a subset of well-measured tracks coming from the beamline to determine the 

374 event's primary vertex. The longitudinal coordinate Zq of the lepton track's point of closest 

375 approach to the beamline must be within 5 cm of the primary vertex to ensure that the 

376 lepton and the jets come from the same hard interaction. We reduce backgrounds from Z 

377 boson decays by vetoing events where the invariant mass of the lepton and a second track 
37B with pt > 10 GeV/c falls in the Z-boson mass window 76 < me-trk < 106 GeV/c^. 

379 We use the b-jet tagging strategy developed in the previous WH search [sl. We require 

380 at least one jet to be 6-tagged with the secondary vertex algorithm, and then we divide our 

381 sample into three exclusive categories of varying purity. Events with two secondary vertex 

382 tagged jets have the highest purity, followed by events with one secondary vertex tagged jet 

383 and one jet probability tagged jet. In the lowest purity events, there is only one secondary 

384 vertex tagged jet. 

385 We further purify the sample with exactly one secondary vertex tagged jet by using 

386 kinematic and angular cuts designed to reject QCD events with fake W signatures. The 

387 kinematics of the QCD contamination vary with the lepton signature they mimic. We 

388 therefore apply a separate veto to each lepton subsample. 



389 One approach we use to reduce QCD is to cut on a varia 



Die correlated with mismea- 



390 surement. The observation of single top quark production [23| demonstrated that missing 

391 transverse energy significance is a useful variable to remove QCD contamination. Miss- 

392 ing transverse energy significance quantifies the likelihood that the measured I^t comes 

393 from jet mismeasurements. is defined as follows: 

^ (Ejets Chs COs'iA<li^rJet)E^T7et + COS^ ( A<^,, w)i^T, w) ' 

395 where Cjes is the jet energy correction factor, Acf^^j^t is the azimuthal angle between the 

396 jet and the I^t direction, is the uncorrected jet Et, unclustered energy is energy 

397 not associated with a jet, Ex^unci is the transverse unclustered energy, and A(j:^^^unci is the 
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Quantity 


Cut 


CEM 


Mt 


> 20 GeV 




> -0.05 • Mt + 3.5 




> 2.5 - 3.125 • A<f>MET,jet2 


CMUP,CMX 


Mt 


> 10 GeV 


ISOTRK 


Mt 


> 10 GeV 



TABLE IV: QCD veto cuts for each lepton category. These cuts are applied to events with exactly 
one identified b-jet. 

398 azimuthal angle between the unclustered energy direction and the I^t direction. The lower 

399 the value of , the more likely it is that the I^t comes from fluctuations in jet energy 

400 measurements. The uncertainty on the calorimeter energy not clustered into one of the jets 

401 is also included. 

402 Another useful approach for rejecting QCD backgrounds is to require that the lepton 

403 momentum and I^t be consistent with the decay of a W boson. However, since only the 

404 transverse component of the neutrino momentum is available via J^t , the W invariant mass 

405 cannot be calculated. Instead, if we ignore the neutrino pz, we can calculate the transverse 

406 mass as follows: 

407 Mt = ^2{pt^^t^pt^^p ■ (3) 

408 We use both Mt and to remove QCD events from our sample. Table [IV] lists the 

409 different QCD veto cuts for each lepton type. The cuts were chosen to have high efficiency 

410 for events with a W boson while rejecting the maximum amount of QCD and minimizing 

411 disagreement between data and MC in the pretag sample. 
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412 IV. BACKGROUNDS 



413 The signature of WH associated production is shared by a number of processes that 

414 can produce the combination iubb. The dominant backgrounds are W^+jets production, tt 

415 production, single top production, and QCD multijet production. Diboson production and 

416 Z+jets production, collectively referred to as "electroweak backgrounds," contribute to the 

417 sample at smaller rates than any of the other backgrounds. Diboson production has a small 
41B contribution because of its small cross section and, in the case of WW, lack of 6-jets at leading 

419 order. Z+jets production has a small contribution because it has a small overlap with our 

420 single lepton final state. Our estimate of the background rates uses a combination of Monte 

421 Carlo techniques and data-driven estimates. Our data-driven estimates use background- 

422 enriched control regions outside of our search region to determine background properties. 

423 We extrapolate the background properties from the control regions to the search region and 

424 assess an uncertainty on the estimates. Our background techniques are common to top cross 

425 section measurements 



21| , single top searches [2j] , and prior WH searches [25| . We provide 

426 an overview of the background estimate below and discuss the details of each background 

427 in the subsections that follow. 

42B We first describe our background estimate for the sample of iujj events without any tag- 

429 ging requirements applied, which we refer to as the pretag sample. This sample is composed 

430 of events from two classes of processes: (1) events containing a high-p^ lepton from a real W 

431 decay and (2) events in which the lepton is from a source other than a W. In the second class 

432 of events, referred to as QCD multijet events, the high-p^- lepton comes either from a jet 

433 that fakes a lepton signature or from a real lepton produced in a heavy-flavor decay. After 

434 the QCD multijet background is subtracted off, what remains are events from a collection 

435 of processes that include the production of a.W boson: primarily W + jets production, top 

436 production, and other electroweak backgrounds. We use a Monte Carlo based technique 

437 to estimate the relative contributions of processes whose rates and topologies are described 
43B well by next-to-leading order (NLO) calculations. These processes include ti, single top and 

439 diboson, and Z + jets production. We estimate their expected contribution N using the 

440 theoretical NLO cross section a, Monte Carlo event detection efficiency corrected to match 

441 the efficiency in the data e, and the integrated luminosity of our dataset Cint- 

N = a-e- Cm (4) 
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443 We subtract the contribution of these processes from the total number of observed events. 

444 After accounting both for the fraction of QCD multijet events and for the top and other 

445 electroweak processes, what remains are the pretag PF+jets events, whose contribution is 

446 estimated as follows: 

w ^W+Jets = Npretag ' (1 " Fqcd) - NewK " NtOP (5) 

448 where Npretag is the observed number of iujj pretag events, Newk is the number of esti- 

449 mated electroweak events, and Ntop is the number of estimated top events. 

450 We estimate the number of tagged W + jets events using the number of pretag W + jet 

451 events and a tag probability. Wc measure the tag probabilities for both light and heavy- 

452 flavor jets in inclusive jet data. The tag probability for heavy-flavor jets is etag, and the tag 

453 probabihty for falsely tagged jets, called "mistags" , is emistag- W + bb, W + cc, and W + cq 

454 production are collectively referred to as W + heavy-flavor processes. All other W + jets 

455 production is referred to as -|- light flavor. We use a 6-tag scale factor to correct the 

456 Monte Carlo tagging efficiency according to the tag efficiency observed in data. We must 

457 estimate the fraction of + jet events that are W + heavy-flavor events Fhf in our sample 

458 in order to use the appropriate tag probabilities. We use -|- 1 jet data to calibrate the 

459 heavy-flavor fraction from the Monte Carlo. We use the ratio of the heavy-flavor fraction in 

460 the data F^p" to the heavy-flavor fraction in the Monte Carlo F^f to calculate a correction 

461 factor K — F^p-fF^p. We apply the correction factor to the number of 1^ -|- heavy jets 

462 estimated with the Monte Carlo. After including this cahbration, the number of l^-|-jets in 

463 the tagged sample is: 

N'wThf = N^tTets ■ {Fhf ■ K) ■ e,„, (6) 

'^65 ^W+LF = ^W+jets ■ {i - Fhf ■ K) ■ ^mistag (7) 

466 The estimation of the rate of these backgrounds are done separately for each jet bin in 

4«7 the data. Below we describe the estimation of the individual pieces in greater detail. 
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Process Theoretical Cross Section 



WW 
WZ 

zz 

Single top s-channel 
Single top t-channel 
ti 

Z + Jets 



TABLE V: Theoretical cross sections 



12.40 ± 0.80 pb 
3.96 ± 0.06 pb 
1.58 ± 0.05 pb 
0.88 ± 0.11 pb 
1.98 ± 0.25 pb 
6.7 ± 0.83 pb 
787.4 ± 85 pb 



26l428l| and uncertainties for the electroweak and top 
backgrounds. Top cross sections assume a mass of rrit = 175GeV/c^. 



A. Top and Electroweak Backgrounds 



469 The normalization of the diboson, Z+_ 



on the theoretical cross sections 
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26 



jets , top-pair, and single-top backgrounds are based 
28[ listed in Table |Vl The estimate from theory is 



471 well-motivated because the cross sections for most of the processes have small theoretical 

472 uncertainties. Z-|-jets is the only process where the large corrections to the leading order 

473 process give large uncertainties to the theoretical cross section. The impact of the large 

474 uncertainty on our sensitivity is marginalized by the small overlap of Z-l-jets with the VT-l-jets 

475 final state. The background contributions are estimated using the theory cross sections, 

476 luminosity, and the Monte Carlo acceptance and 6-tagging efficiency. The Monte Carlo 

477 acceptance is corrected for lepton identification, trigger efficiencies, and the z vertex cut. 

478 We also use a b-tagging scale factor to correct for the difference in tagging efficiency in 

479 Monte Carlo compared to data. 



480 B. QCD Multijet 

481 QCD multijet events can fake a W signature when a jet fakes a lepton and overall mis- 

482 measurement leads to fake Ifx ■ Since these events do not have real W bosons in them, we 

483 also use the term non-W to refer to QCD multijet events. It is difficult to identify the precise 

484 sources of mismeasurement and handle them appropriately in a detector simulation. The 
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485 difficulty is increased by the large number of processes that contribute to the composition 

486 of the QCD background at unknown relative rates. Each lepton category is susceptible to 

487 different kinds of fakes. We use different QCD models for central-lepton triggered events 

488 and isolated track events. 

489 We model central-lepton triggered QCD events using events where a jet ffied the electron 

490 trigger, passed the electron kinematic cuts, but failed exactly two of the calorimeter or 

491 tracking quality cuts. Events that fail these cuts will have the kinematic properties of W 

492 events, including isolation, but the sample will be enriched in fakes. This is the same model 



493 used in the CDF observation of single top 



23| . As noted in that paper, these fake events 



494 have the remarkable property that they model both electron and muon fakes. 

495 We model QCD events that fake an isolated track by using events recorded on the I^t + 

496 2 Jets trigger. We use events with muon candidates that are not calorimeter isolated and 

497 are within the isolated track acceptance {\ri\ < 1.2). Calorimeter isolation is defined as the 

498 fraction of the lepton energy in a cone of AR = 0.4 surrounding the lepton. Non-isolated 

499 leptons are unlikely to come from the decay of an on-shell W, and thus are enriched in fakes. 

500 We estimate the amount of QCD background in each sample by fitting the I^t spectrum 

501 in data. The fit includes the control region i^^ < 20 GeV, which is enriched in QCD fakes. 

502 Figure |2] shows the I^t fit for isolated track pretag events. The fit has one component with 

503 fixed normalization and two templates whose normalizations can vary. The fixed component 

504 is a combination of top and electroweak processes whose normalizations are described in 

505 Section IIV A[ We let the W + jets template vary along with the QCD template because 

506 there is a large uncertainty on the VT+jets cross section. The QCD template has a I^t 

507 spectrum that peaks near low I^t , and its normalization is driven by the low i^j- bins. The 

508 normalization of the ly+jets template is driven by the high I^t region. The fit determines 

509 the relative amounts of QCD and VT+jets in the full I^t sample, and we use these fit results 

510 to determine the QCD fraction in the search region {I^t > 20 GeV). For isolated track events 

511 with two jets and no b-tag requirement, we estimate a 19% QCD fraction in the signal region, 

512 as shown in Fig. |2l The pretag QCD fractions for the other lepton types are less than the 

513 isolated track fractions. Pretag CEM electrons events have 10% QCD fraction, and both 

514 CMUP and CMX muon events have a 3% QCD fraction. While isolated tracks have a larger 

515 amount of QCD events than the other lepton types, the vast majority of the isolated track 

516 events (81%) still contain W bosons. We use the QCD fractions for each lepton type and 
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517 tag category in the calculations for the background summaries in Tables IVIIII through IXIII] 
51B We estimate the uncertainty of the QCD normalization by studying the change in the 

519 QCD fraction due to changes in the QCD model. For tight lepton events we use an alternate 

520 QCD model based on leptons that fail our isolation requirements. We find a 40% uncertainty 

521 to the QCD normalization that covers the effect of using this alternative model. We use the 

522 same uncertainty estimate for both tight leptons and isolated tracks. 



CDF RUN II 2.7 fb 




FIG. 2: Fit of the pretag isolated track Ffx control region that is used to determine the QCD 
fraction of isolated track events. The arrow illustrates the FfT cut. We estimate a QCD fraction of 
19% for the region with Ffx > 20 GeV. There is some disagreement between the data and our model 
in the Iow-I^t control region, and also around 50-55 GeV. The figure shows just one QCD model. 
The difference between this nominal model and are alternate covers the modelling difference shown 
here. We use the difference between the two models as our systematic uncertainty. 



523 C. ly+Heavy-Flavor 

524 The number of W + heavy flavor events is a fraction the number of W + light flavor 

525 events, as described by Fhf in Equations [6] and [71 The fraction of l^+heavy-flavor events has 



526 been studied extensively and is modeled in the ALP GEN Monte Carlo Generator 



30|. We 



527 calibrate the ALPGEN Version 2 W + jets Monte Carlo heavy-flavor fraction to match the 
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observed heavy flavor fraction in the 1^ + 1 jet control region. We use the same cahbration 



231]. The cahbration uses template 



529 of the heavy-flavor fraction as the single top observation 

530 flts of flavor-separating variables in 6-tagged + 1 jet data to measure the heavy flavor 

531 fraction. The calibration measures the calibration factor as deflned in equation [HI to be 

532 K =IA± 0.4. 

533 We can estimate the amount of + heavy flavor events in our signal region by calculating 

534 the efficiency for these events to pass our tag requirements etag- The efficiency etag is 

jets 

etaa = 1 - n(l-Ptag), (8) 

i 

536 where p\^g is the probability for jet i in the event to have a 6-tag. The probability for a h- 

537 tagged Monte Carlo jet originating from a 6 or c quark to have a 6-tag in the data is the 6-tag 

538 scale factor. The 6-tag scale factor is the ratio of data to Monte Carlo 6-tag efficiencies. It is 

539 estimated to be 0.95 ± 0.04 for secondary vertex tags [3] and 0.85 ± 0.07 for jet probability 

540 tags [22]. In the case where there are additional light-flavor jets produced in the W + heavy 

541 flavor events, there is a small chance for those light-flavor jets to be incorrectly tagged as 

542 6-jets. We account for this in Equation [H] by giving these just a small probability to be 

543 incorrectly tagged. We call the probability to be incorrectly tagged the mistag probability. 

544 It is discussed in detail in Section IIVDI 

545 Table |VT] shows the corrected heavy-flavor fractions for our W + heavy-flavor samples di- 

546 vided according to the heavy-flavor process and number of reconstructed jets. It is necessary 

547 to divide the samples by heavy-flavor process because h- and c-jets have different tagging 
54B efficiencies. Table IVIII shows the corrected per-event tagging efficiencies. We calculate the 

549 W + heavy-flavor normalizations using Eq. [6] and the fractions and efficiencies from the 

550 tables. 

551 The two sources of uncertainties for the W + heavy-flavor backgrounds are the 6-tag 

552 scale factor uncertainty and the heavy flavor fraction uncertainty. We accommodate the 

553 6-tag scale factor uncertainty by shifting the scale factor by ±lo", propagating the change 

554 through our background calculation, and using difference between the shifted and nominal 

555 calculation as our error. We add this error in quadrature with the heavy-flavor fraction 

556 uncertainty and use the total error as a constraint on the background in our likelihood flt. 
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Corrected Heavy Flavor (HF) fraction (%) 
of inclusive + jet events by jet multiplicity 



Process Number of Jets Fraction of Events by Jet Multiplicity 

matched to HF + 2 jets 1^ + 3 jets W + A jets 1^ + 5 jets 
Wbb (16) 2.2 ± 0.88 3.5 ± 1.4 4.63 ± 1.8 5.5 ± 2.2 

Wbb (26) 1.32 ± 0.52 2.6 ± 1.0 4.17 ± 1.7 6.0 ± 2.4 

Wcc (Ic) 11 ± 4.4 14 ± 5.6 15.18 ± 6.1 15.8 ± 6.3 

W^cc (2c) 2.1 ± 0.84 4.7 ± 1.9 7.69 ± 3.1 10.9 ± 4.4 

TABLE VI: The corrected fraction of inclusive W + jet events that contain heavy-flavor. The 
fractions are divided into separate categories according to the Monte Carlo flavor information for 
jets in the event and the number of reconstructed heavy-flavor jets. For example, Wbb (16) events 
have two 6-quarks at the generator level, but only one 6-quark matched to a reconstructed jet. 
The fractions from ALPGEN Monte Carlo have been scaled by the data-derived calibration factor 
of 1.4 ±0.4. 



557 D. Mistagged Jets 

558 W + light flavor events with a fake 6-tag migrate into our signal region. Our estimate 

559 of the number of falsely tagged ly+light flavor events is based on the pretag number of W 

560 + light flavor events and the sample mistag probability Cmistag in equation [3 The sample 

561 mistag probability is based on the per-jet mistag probability. For each event in our W + 

562 light flavor Monte Carlo samples, we apply the per-jet mistag probability to each jet and 

563 combine the probabilities to get an event mistag probability. We combine the event mistag 

564 rates to get emistag- 

565 We estimate the per-jet mistag probability for each of our two tagging algorithms using 

566 a data sample of generic jets with at least two well-measured silicon tracks. The decay 

567 length is defined as the distance between the secondary vertex and the primary vertex in the 

568 plane perpendicular to the beam direction. This decay length is signed based on whether 

569 the tracks are consistent with the decay of a particle that was moving away from (positive 

570 sign) or towards (negative sign) the primary vertex. False tags are equally likely to have 

571 positive or negative decay lengths to first order. The symmetry allows calibration of the 
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Corrected Per-event 5-tag efficiencies 



One SECVTX Tag Efficiency 



Jet Multiplicity 


2 jets 


3 jets 






Event Eff (16) (%) 


23.10 


24.68 


25.02 


27.14 


Event Eff (26) (%) 


30.09 


30.34 


30.35 


29.71 


Event Eff (Ic) (%) 


7.02 


7.69 


8.68 


10.24 


Event Eff (2c) (%) 


9.46 


10.46 


11.24 


12.12 


Two SECVTX Tag 






Jet Multiplicity 


2 jets 


3 jets 


4 jets 




Event Eff (16) (%) 


0.30 


0.78 


1.34 


1 76 


Event Eff (26) (%) 


8.76 


9.68 


10.18 


11.14 


Event Eff (Ic) (%) 


0.04 


0.12 


0.24 


0.40 


Event Eff (2c) (%) 


0.38 


0.55 


0.88 


0.91 


One SECVTX TAG + 


ne JETPROB Ta; 


5 Efficiency 


Jet Multiplicity 


2 jets 


3 jets 


4 jets 


5 jets 


Event Eff (16) (%) 


0.79 


1.75 


2.57 


3.74 


Event Eff (26) (%) 


6.95 


7.78 


8.86 


9.77 


Event Eff (Ic) (%) 


0.20 


0.47 


0.78 


1.24 


Event Eff (2c) (%) 


1.19 


1.59 


2.14 


2.43 



TABLE VII: The corrected per-event tagging efficiencies for events with heavy-flavor content. The 
event efficiencies are divided into separate categories depending on the Monte Carlo truth flavor 
information for jets in the event: 16 events have one jet matched to 6-quark, 26 events have two 
jets matched to a 6-quark, Ic events have one jet matched to a c-quark, and 2c events have two 
jets matched to a c-quark. 



572 false tag probability using negative tags. There is a slightly greater chance for a false tag 

573 to have a positive decay length due to material interaction, and our estimate accounts for 

574 this asymmetry. The false tag probability for secvtx is parameterized in bins of ?7, number 

575 of vertices, jet E^, track multiplicity, and the scalar sum of the total event Et 21|. We 

576 parameterize jet probability mistaging in jet r/, z position of primary vertex, jet E^, track 
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577 multiplicity, and scalar sum of the total event Et- 

57B We estimate the uncertainties on the per-jet mistag probability by using negatively tagged 

579 jets in the data. The uncertainty estimates check for consistency between the number of 

580 expected and observed negative tags. The uncertainties are accounted for in the analysis by 

581 fluctuating the per-jet tag probabilities by ±la, and propagating the change through the 

582 background estimate. 

583 E. Summary of Background Estimate 

584 Tables IVIIII through IXIII] summarize our background estimate for our dataset of 2.7 fb~^. 

585 Figures [3] through |5] present the information from the tables as plots. The plots show the 

586 background estimate compared to data. The largest errors on the background estimate come 

587 from the large uncertainty on the heavy flavor fraction used to calculate W + charm and 

588 W + bottom. We add these large uncertainties linearly because they come from the same 

589 source. The 6-tagging scale factor uncertainty is also correlated across all backgrounds and 

590 added linearly. In general, the background estimate agrees with the data within uncertainties 

591 for each jet multiplicity. The agreement of the background estimate with the data in the 

592 high-jet-multiplicity bins gives us confidence that our estimate is correct in our two-jet search 

593 region. 
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CDF Run II 2.7 fb-^ 
Tight Lepton Background Prediction and Event Yields 
Events with Exactly One Secvtx Tag 



Pf( )COSS 


2j('(.s 


3j('is 


lj(Ms 


jets 


All P*rpt;fi0" (HfiTifliflflt.p^^ 


38729 


6380 


1677 


386 


WW 


40.6 ± 4.2 


11.9 lb 1.2 


2 92 ± 25 


71 ± 06 


wz 


18 86 ± n 94 


3 43 ± 23 


93 ± 06 


2 ± 02 


zz 


fl 48 ± 04 


19 ± 06 1 


081 ± 007 


023 ± 002 


Ton Pair 


102 ± 14 


193 ± 26 


183 ± 26 


59 4 ± 8 8 


Si 1 n crl *3 T'r^r^ q_ (; n q n n f^l 






1 47 _|_ n 1 c: 




Sirip"lp T'oD t-C^liannpl 


42 53 ±44 


9 24 lb 94 


1 62 lb 17 


22 lb 02 


Z+Jets 


28.72 ± 3.4 


8.65 lb 0.96 


2.73 lb 0.29 


0.53 lb 0.06 


W+bottom 


365.6 ± 140 


91.0 ± 35 


19.4 lb 8 


3.97 ± 1.7 


W+charm 


364.6 ± 140 


81.2 lb 31 


17.3 ± 7 


3.64 ± 1.6 


Mistags 


319 ± 42 


83.8 ± 13 


18.8 ± 5.07 


3.82 lb 1.5 


Non-W 


107 ± 43 


40.2 lb 17 


17.3 lb 14 


4.48 lb 4.4 


Total Prediction 


1408 ± 287 


530 lb 75 


266 ± 34 


77 lb 11 


ObscM'ved 


1404 


■186 


281 


81 



TABLE VIII: Background summary table for events with a central lepton and exactly one secondary 
vertex tag. The heavy-flavor fraction F^p is the source of the large correlated uncertainty for 
W^+bottom and W^+charm. The other large source of correlated uncertainty is the 6-tagging scale 
factor. 
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CDF Run II 2.7 fb-^ 
Isolated Track Background Prediction and Event Yields 
Events with Exactly One Secvtx Tag 



Pt'OCX'SS 


2jets 


3j(Ms 






All Prpf.fiP" (HfiTifliflflt.p^ 


4253 


1380 


427 


117 


WW 


6.4 ± 0.65 


2.83 ± 0.25 


0.75 ± 0.07 


0.23 lb 0.02 


wz 


2.41 ± 0.16 


0.92 ± 0.06 


0.19 ± 0.01 


0.063 zb 0.005 


TIL 


0.127 ± 0.009 0.052 ± 0.004 0.007 ± 0.001 0.006 ± 0.001 


Ton Pair 


28.0 ± 3.8 


58.3 ± 8.0 


53.4 ± 7.6 


16.8 lb 2.5 


^in cflf» T'r^Tv G-l; n an ■n^il 


6.08 ± 0.58 


1.91 ± 0.19 


0.43 ± 0.04 


0.08 ± 0.01 


SiTiP'lp Ton t-(~!bariTipl 


10.1 ± 1.1 


2.32 ± 0.24 


0.41 ± 0.05 


0.07 lb 0.01 


Z+Jets 


9.05 ± 1.1 


3.35 ± 0.36 


0.74 ± 0.077 


0.16 zb 0.02 


W+bottom 


39.9 ± 16 


18.4 ± 7.3 


5.35 ± 2.3 


1.91 ± 0.79 


W+charm 


36.7 ± 15 


16.2 ± 6.5 


4.66 ± 2.0 


1.53 ± 0.64 


Mistags 


43.2 ± 8.2 


17.7 ± 4.0 


4.81 ± 1.7 


1.82 zb 0.64 


Non-W 


37.6 ± 15 


22.2 ± 8.9 


5.26 ± 4.2 


2.13 zb 1.7 


Total Prediction 


220 lb 35 


144 ± 19 


76 lb 10 


25 zb 3.4 


01)sorA"od 


208 


150 


78 


31 



TABLE IX: Background summary table for events with an isolated track and exactly one secondary 
vertex tag . The heavy-flavor fraction F^f is the source of the large correlated uncertainty for 
W^+bottom and VK+charm. The other large source of correlated uncertainty is the 6-tagging scale 
factor. 
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CDF Run II 2.7 fb-^ 

Tight Lepton Background Prediction and Event Yields 
Events with One Secvtx Tag and One Jet Prob Tag 



Pt'OCX'SS 


-j^ ' 


3 jots 


lj(Ms 


5 jets 


All Prpf.fiP" (HfiTifliflflt.p^ 


44723 


7573 


1677 


386 


WW 


1 24 ± n 58 


n 85 + 81 


0.4 lb 0.13 


165 ± 047 


wz 


2 51 -t n 48 


78 ± 16 


18 ± 04 


052 ± 013 


TIL 




053 + 009 


021 + 004 


f) one _|_ f) nni 


Ton Pair 


20.4 ± 4.2 


63 9 ± 13 


79.3 ± 16 


29 9 ± 6 1 


Sin gtIo TWtv g-I; n an ■n^il 


« qq _l_ 1 1 




f) C7 _|_ n 1 


n 1 QO _|_ f) 094 


SiTio'lp Tod t,-(~]bariTipl 


2.1 ± 0.64 


1 67 lb 86 


46 lb 09 


076 lb 015 


Z+Jets 


1.81 ± 0.54 


1.17 lb 0.35 


0.34 lb 0.12 


0.1 lb 0.03 


W+bottom 


49.1 ± 20 


17.1 lb 7.2 


4.89 lb 2.1 


1.28 ± 0.59 


W+charm 


18.0 ± 8.3 


7.89 ± 3.7 


2.57 lb 1.2 


0.67 ± 0.34 


Mistags 


5.84 ± 6.0 


3.01 lb 3.4 


0.1 lb 1.1 


0.29 ± 0.37 


Non-W 


11.1 ± 5.38 


6.57 lb 3.5 


8.88 lb 3.4 


1.51 lb 2.1 


Total Prediction 


119 lb 30 


105 lb 19 


93 lb 17 


34 lb 7 


01)sorA"od 


124 


109 


101 


36 



TABLE X: Background summary table for events with a central lepton and two tags: one sec- 
ondary vertex tag and one jet probability tag. The heavy-flavor fraction Fhf is the source of the 
large correlated uncertainty for W^+bottom and VF+charm. The other large source of correlated 
uncertainty is the 6-tagging scale factor. 
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CDF Run II 2.7 fb-^ 

Isolated Track Background Prediction and Event Yields 



Events with One Secvtx Tag, One 


Jet Prob Tag 






2jets 




lj(Ms 




All Prp'i'PiO" (^fmHiHflf P<3 


5149 


1623 


487 


124 


WW 

V V V V 


0.2 ± 0.09 


0.24 ± 0.09 


0.1 ± 0.03 


0.03 lb 0.01 


wz 

V V ZJ 


0.51 ± 0.09 


0.2 ± 0.04 


0.048 ± 0.01 


0.013 zb 0.004 


7.7. 


0.032 ± 0.006 0.021 ± 0.005 0.007 ± 0.001 0.002 ± 0.001 


Ton Pair 


6.44 ± 1.3 


20.0 ± 4.2 


24.6 ± 4.9 


8.98 ± 1.8 


oniric ±UJJ o~v^liclliilci 


1.93 ± 0.31 


0.74 ± 0.13 


0.19 ± 0.03 


0.043 ± 0.009 


SiTiP'lp ToTi f-Clnarmpl 


0.53 ± 0.16 


0.5 ± 0.11 


0.12 ± 0.03 


0.028 lb 0.005 


Z+Jets 


0.61 ± 0.2 


0.41 ± 0.13 


0.13 ± 0.04 


0.039 lb 0.013 


W+bottom 


6.0 ± 2.7 


3.4 ± 1.6 


1.37 ± 0.67 


0.59 ± 0.26 


W+charm 


2.14 ± 1.07 


1.64 ± 0.86 


0.77 ± 0.41 


0.34 ± 0.17 


Mistags 


0.8 ± 1.18 


0.61 ± 0.84 


0.27 ± 0.31 


0.13 ± 0.17 


Non-W 


1.97 ± 0.79 


1.38 ± 0.55 


0.99 ± 0.79 


0.37 lb 0.5 


Total Prediction 


21 ± 4 


29 lb 5 


29 lb 5 


11 lb 2 


01)sorA"od 


21 


30 


32 


12 



TABLE XI: Background summary table for events with an isolated track and two tags: one sec- 
ondary vertex tag and one jet probability tag. The heavy-flavor fraction Fhf is the source of the 
large correlated uncertainty for W^+bottom and VF+charm. The other large source of correlated 
uncertainty is the 6-tagging scale factor. 
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CDF Run II 2.7 fb-^ 
Tight Lepton Background Prediction and Event Yields 





Events with Two Secvtx 


Tags 




Process 


2 jets 


3i(Ms 


Ijcis 


5j(-ls 


All Pretag Candidates 


4472 S 


7573 


1677 


386 


WW 


n Q _|_ n n« 


n oq + 05 


0.17 ± 0.03 


0.08 ± 0.01 


wz 


9 Qo + n S7 

kJ.kJ^ —1— \J-kJ I 


n q4 _)_ f) 1 1 


0.19 ± 0.02 


0.04 ± 0.01 


zz 


\j.± _i_ \j.\j± 


n (-170 _|_ n nno 


0.019 ± 0.002 


0.005 ± 0.001 


Top Pair 




7fi 8 + 1 2 


101 ± 16 


36.1 ± 5.9 


Single Top s-Channel 


q cjc _|_ 1 2 


Q 25 + n 41 


0.72 ± 0.09 


0.15 ± 0.02 


Single Top t-Channel 


2 15 + 03 


1 9 ± 26 


0.53 ± 0.07 


0.1 ± 0.01 


Z+Jets 


1.42 ± 0.2 


0.95 ± 0.13 


0.26 ± 0.04 


0.085 ± 0.013 


W+bottom 


55.0 ± 22 


18.1 ± 7.4 


4.88 ± 2.0 


1.24 ± 0.55 


W+charm 


4.87 ± 2.0 


2.35 ± 1 


0.94 ± 0.4 


0.25 ± 0.12 


Mistags 


1.38 ± 0.39 


0.93 ± 0.3 


0.34 ± 0.12 


0.11 ± 0.05 


Non-W 


8.96 ± 4.0 


5.02 ± 2.0 


0.74 ± 1.6 


0.23 ± 1.5 


Total Prediction 


113 ± 25 


111 ± 16 


110 lb 17 


38 lb 6 


ObscrA'od 


114 


132 


104 


42 



TABLE XII: Background summary table for events with a central lepton and two secondary vertex 
tags. The heavy-flavor fraction Fhf is the source of the large correlated uncertainty for 1^+bottom 
and l^+charm. The other large source of correlated uncertainty is the 6-tagging scale factor. 
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CDF Run II 2.7 fb-^ 
Isolated Track Background Prediction and Event Yields 
Events with Two Secvtx Tags 



Pt'OCX'SS 


-j^ ' 


3 




s 




5 jc'ls 


All Prpf.fiP" (nfiTifliflflt,P9 


5149 


1623 


487 


124 


WW 




U. io 


± 


n n9 
u.uz 


067 ± 012 


019 ± 003 


wz 


65 ± n (17 


n 9/1 


It 


u.uo 


029 ± 003 


01 ± 001 


TIL 


045 + 005 


U.Uzo 


1 

± 


U.UUo 


01 + 001 


002 + 

\j .\j\J £j _l_ \J 


Ton Pair 


7.75 ± 1.2 


09 7 


■ _i_ 


Q 7 
o. < 


31.5 ± 5.1 


11.5 ± 1.9 


^ 1 n cf 1 f» T'r^T^ G- 1! n Q n n fil 




0.91 


± 


0.12 


n 91 _|_ 03 


045 + 006 


SiTiP'lp Trfn f-linj^Tinpl 


58 ± n (18 


0.57 




0.08 


18 ± 02 


035 ± 005 


Z+Jets 


0.51 ± 0.07 


0.32 


± 


0.05 


0.093 ± 0.014 


0.025 ± 0.004 


W+bottom 


7.51 ± 3.3 


3.59 




1.63 


1.41 ± 0.66 


0.53 ± 0.23 


W+charm 


0.68 ± 0.3 


0.56 




0.26 


0.26 ± 0.13 


0.18 ± 0.05 


Mistags 


0.27 ± 0.13 


0.2 




0.1 


0.089 ± 0.05 


0.052 ± 0.026 


Non-W 


1.78 ± 0.71 


1.89 




0.76 


6.53 ± 5.2 


2.65 ± 2.1 


Total Prediction 


22 ± 4 


31 


± 


4 


40 ± 7 


15 lb 3 


Observed 


24 




31 




37 


15 



TABLE XIII: Background summary table for events with an isolated track and two secondary 
vertex tags. The heavy-flavor fraction F^p is the source of the large correlated uncertainty for 
VF+bottom and 14^+charm. The other large source of correlated uncertainty is the 5-tagging scale 
factor. 
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CDF Run II 2.7 fb"^ 
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FIG. 3: Number of expected and observed background events for events with exactly one secvtx 
tag, shown as a function of jet multiphcity. The plots show tight leptons (top) and isolated tracks 
(bottom). The hatched regions indicate the total uncertainty. 
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CDF Run II 2.7 fb"^ 
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CDF Run II 2.7 fb' 
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FIG. 4: Number of expected and observed background events for events with one secvtx tag 
and one jetprob tag, shown as a function of jet multiplicity. The plots show tight leptons (top) 
and isolated tracks (bottom). The hatched regions indicate the total uncertainty. 
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CDF Run II 2.7 fb"^ 
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CDF Run II 2.7 fb' 
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FIG. 5: Number of expected and observed background events for events with two secvtx tags, 
shown as a function of jet multiphcity. The plots show tight leptons (top) and isolated tracks 
(bottom). The hatched regions indicate the total uncertainty. 
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594 V. HIGGS BOSON SIGNAL ACCEPTANCE 



595 We simulated the WH signal kinematics using the pythia Monte Carlo program [31 1. 

596 We generated signal Monte Carlo samples for Higgs masses between 100 and 150GeV/c^. 

597 The number of expected WH — )■ iubb events, A^, is given by: 

598 N = €■ J Cdt-a{pp-^WH)-B{H ^bb), (9) 

599 where e, / Cdt, a{pp — )■ WH), and B{H — )• bb) are the event detection efficiency, integrated 

600 luminosity, production cross section, and branching ratio, respectively. The production cross 

601 section and branching ratio are calculated to next-to-leading order (NLO) precision j4|. 

602 The total event detection efficiency is composed of several efficiencies: the primary vertex 

603 reconstruction efficiency, the trigger efficiency, the lepton identification efficiency, the b- 

604 tagging efficiency, and the event selection efficiency |5|. Each efficiency is calibrated to 

605 match observations. 

606 We parametrize the Ifx trigger turn-on as a function of ^J?''*®^ , which is Efx corrected for 

607 the primary vertex position but not muons or jet energy scale corrections. We use Ip^^^^^ 

608 because it is close to the Bfr calculation used by the trigger and is modeled better in the 

609 Monte Carlo than which is calculated assuming = 0. The measurement of the jets 

610 can influence the measurement of the Ifx ■ We require that the jets in the event are above 

611 the trigger threshold {Et > 25 GeV) and well separated {AR > 1.0), which reduces the 

612 impact of the jets on the J^t ■ We measured the turn-on curve using events recorded with 

613 the CMUP trigger, which is independent from the I^t + 2 jets trigger. We selected events 

614 passing our jet requirements, and measured their efficiency to pass the J^t + 2 jets trigger 

615 as a function of ^Jf''**^^ . Figure E] shows the measured I^t + 2 jets trigger turn-on. We use 

616 the parmeterized turn-on curve to weight each Monte Carlo event according to its efficiency 

617 to pass the trigger. 

618 The expected number of signal events is estimated by equation |9] at each Higgs boson 

619 mass point. Table IXTVl shows the number of expected WH events for Mh = 120 GeV/c^ in 

620 2.7 fb-^ 

621 The total systematic uncertainty on the acceptance comes from several sources, including 

622 the jet energy scale, initial and final state radiation, lepton identification, trigger efficiencies, 

623 and 6-tagging scale factor. The largest uncertainties come from the 6-tagging scale factor 

624 uncertainty and isolated track identification uncertainty. 
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10 20 30 40 50 60 70 80 90 100 
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FIG. 6: Bfx plus jets trigger turn-on curve parameterized as a function of vertex Ffx ■ The plot 
shows the turn-on curve measured in 2.7 fb~^ of CDF data. 



625 We assign a 2% uncertainty to the CEM, CMUP, and CMX lepton identification efficiency, 

626 and an 8% uncertainty to isolated track identification. The identification uncertainties are 

627 based on studies comparing Z boson events in data and Monte Carlo. 

628 The high pt lepton triggers have a 1% uncertainty on their efficiencies. We measure the 

629 trigger efficiency uncertainty by using backup trigger paths or Z boson events. We measure 

630 a 3% uncertainty for events collected on the Efr + 2 jets trigger by examining the variations 

631 in the Efx turn-on curve in sub-samples with kinematics different from the average sample. 

632 We use the variation in the Ifx turn-on to calculate a variation in signal acceptance, and we 

633 use the mean variation in signal acceptance as our uncertainty. 

634 We estimate the impact of changes in initial and final state radiation by halving and 

635 doubling the parameters related to ISR and FSR in the Monte Carlo event generation 32|. 

636 The difference from the nominal acceptance is taken as the systematic uncertainty. 

637 The uncertainty in the incoming partons' energies relies on the the parton distribution 

638 function (PDF) fits. A NLO version of the PDFs, CTEQ6M, provides a 90% confidence 

639 interval of each eigenvector The nominal PDF value is reweighted to the 90% confidence 

640 level value, and the corresponding reweighted acceptance is computed. The differences 

641 between the nominal and the reweighted acceptances are added in quadrature, and the total 

642 is assigned as the systematic uncertainty Q]- 
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CDF Run II 2.7 fb-^ 
Number of Expected WH {Mh = 120 GeV/c^) Events 
Lepton Type Expected Number of WH events 



Exactly One Secvtx Tag 



CEM 


1.58 ± 0.08 


CMUP 


0.91 ± 0.05 


CMX 


0.44 ± 0.02 


ISOTRK 


0.72 ± 0.07 


Total 


3.65 ± 0.22 


Two Secvtx Tags 


CEM 


0.66 ± 0.07 


CMUP 


0.37 ± 0.04 


CMX 


0.17 ± 0.02 


ISOTRK 


0.36 ± 0.05 


Total 


1.56 ± 0.18 


One Secvtx Tag and One Jet Probability Tag 


CEM 


0.48 ± 0.05 


CMUP 


0.26 ± 0.03 


CMX 


0.13 ± 0.01 


ISOTRK 


0.23 ± 0.03 


Total 


1.10 ± 0.12 



TABLE XIV: Expected number of WH events at a M(H)=120, shown separately for different tag 
categories and lepton types. The lepton types are categorized based on the sub-detector regions. 

643 The uncertainty due to the jet energy scale uncertainty (JES) [l^ is calculated by shift- 

644 ing jet energies in WH Monte Carlo samples by ±lcr. The deviation from the nominal 

645 acceptance is taken as the systematic uncertainty. 

646 The systematic uncertainty on the 6-tagging efficiency is based on the scale factor un- 

647 certainty discussed in Sec. llVCi The total systematic uncertainties for various 6-tagging 

648 options and lepton categories are summarized in Table IXVI 
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Source 



Uncertainty (%) 



Two Secvtx Tags One Secvtx One JetProb tag Exactly One Secvtx 



Trigger Lepton (Isotrk) ID 


~z/o [o.oo/o) 


~z/o (o.oo/oj 


no/ fd QCO/ N 

~z/o [o.oo/o) 


Lepton (MET+Jets) Trigger 


<1% (3%) 


<1% (3%) 


<1% (3%) 


ISR/FSR 


5.2% 


4.0% 


2.9% 


PDF 


2.1% 


1.5% 


2.3% 


JES 


2.5% 


2.8% 


1.2% 


b-tagging 


8.4% 


9.1% 


3.5% 


Total (Isotrk) 


10.6% (13.8%) 


10.5% (14.0%) 


5.6% (10.1%) 



TABLE XV: Systematic uncertainty on the WH acceptance. "ST+ST" refers to double secondary 
vertex tagged events while "ST+ JP" refers to secondary vertex plus jet probability tagged events. 
Effects of limited Monte Carlo statistics are included in these values. 

649 VI. NEURAL NETWORK DISCRIMINANT 

650 To further improve the signal to background discrimination after event selection, we 

651 employ an artificial neural network (NN). Neural networks offer an advantage over a single- 

652 variable discriminants because they combine information from several kinematic variables. 



3ackgrounds. We 
sl . The following 



653 Our neural network is trained to distinguish W^+Higgs boson events from 

654 employ the same neural network that was used to obtain the 1.9 fb~^ result 

655 section reviews its main features. 

656 Our neural network configuration has 6 input variables, 11 hidden nodes, and 1 output 

657 node. The input variables were selected by an iterative network optimization procedure 
65B from a list of 76 possible variables. The optimization procedure identified the most sensitive 

659 one-variable NN, then looped over all remaining variables and found the most sensitive two- 

660 variable NN. The process continued until adding a new variable does not improve sensitivity 

661 by more than 0.5 percent. The 6 inputs are: 

662 Mjj+: The dijet mass plus is the invariant mass calculated from the two reconstructed jets. 

663 If there are additional loose jets present, where loose jets have Et > 12 GeV, \ri\ < 2.4 

664 and have a centroid within AR < 0.9 of one of the leading jets, then the loose jet that 

665 is closest to one of the two jets is included in this invariant mass calculation. 
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666 J2 Er^Loose Jets): This variable is the scalar sum of the loose jet transverse energies. 

667 pt Imbalance: This variable expresses the difference between I^t and the scalar sum of 

668 the transverse momenta of the lepton and the jets. Specifically, it is calculated as 
Prijeh) + PT{jet2) + Prilep)- ■ 

670 M/™": This is the invariant mass of the lepton, I^j- , and one of the two jets, where the jet 

671 is chosen to give the minimum invariant mass. For this quantity, the Pz component of 

672 the neutrino is ignored. 

673 AR{\epton-h'max)'- This is the AR separation between the lepton and the neutrino. We 

674 calculate the Pz of the neutrino by constraining the lepton and the I^t to the W mass 

675 (80.42 GeV/c^). The constraint produces a quadratic equation for pz and we choose 

676 the larger solution. 

677 Pt{W + H): This is the total transverse momentum of the W plus two jets system, PT{lep+ 

678 + jeti + jet2). 

679 The strongest discriminating variable in the neural network is the dijet mass plus. 

680 We train our neural network with l^+jets, tt, single top, and WH signal Monte Carlo. 

681 We do not use QCD events to train our neural network. We use the same topology and 

682 input variables to train separate neural networks for each Higgs signal Monte Carlo sample. 

683 The samples range from M{H) = 100 to 150 GeV/c^ in 5 GeV increments. At each Higgs 

684 mass, we use the same neural network for tight lepton and isolated track events. 

685 Figures through |9] show the six neural network input variables for isolated track events 

686 in the pretag control region. The plots show that our background model describes the data 

687 reasonably for all the neural network input variables. The modeling is not ideal in regions 

688 that have a large amount of QCD, such as the region around ARmax{MET,1) = 2.5 in 

689 Figure [9] and the region around M^™" = 50 in Figure [HI Figures [10] through [12] show 

690 that these differences are less significant after removing some of the QCD contamination 

691 with 6-tagging. The hashed region in Figures [ID] through [12] indicates uncertainty on the 

692 background estimate. Taking into account the uncertainty on the background estimate, this 

693 modeling is reasonable for the isolated track neural network input variables. 

694 We studied the impact of QCD shape modeling in the tight lepton sample. We did 

695 not expect the QCD shape to have a large impact on the sensitivity because the neural 
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696 network was not trained with QCD events. We found that the large QCD normahzation 

697 uncertainty (40%) accounted for the small variations that arose from using an alternative 

698 QCD model with different kinematics. Based on the tight lepton studies, we assume that 

699 the impact of QCD shape modeling on isolated track sample is also small compared to the 

700 QCD normalization uncertainty. This is not an aggressive assumption since the isolated 

701 track sample only accounts for 20% of the total sensitivity. 

702 The tight lepton categories also show good agreement with the previous publication 
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FIG. 7: Neural network input distributions for isolated track + 2 jet events in the pretag control 
region. The distributions shown are Mjj^ (left) and ^i?T(Loose Jets) (right). The differences in 
shape are attributable to QCD and are less significant in our higher-purity search regions. 



703 VII. LIMIT ON HIGGS BOSON PRODUCTION RATE 

704 We search for an excess of Higgs signal events in our neural network output distributions 

705 using a binned likelihood technique. Figures [T3] through [T5] show the neural network output 

706 distributions for events in different lepton and tag categories. We use the same likelihood 

707 expression and maximization technique as the prior CDF result and described in 34 1- We 

708 maximize the likelihood, fitting for a combination of Higgs signal plus backgrounds. We find 

709 no evidence for a Higgs boson signal in our sample, and so we set 95% confidence level upper 

710 limits on the WH cross section times branching ratio: a{pp — )■ W^H) ■ B{H — t- hb). We 

711 compare our observed limits to our expected sensitivity by creating pseudo-experiments with 
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FIG. 8: Neural network input distributions for isolated track W + 2 jet events in the pretag 
control region. The distributions shown are MJ^J^ (left) and Pt Imbalance (right). The differences 
in shape are attributable to QCD and are less significant in our higher-purity search regions. 
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FIG. 9: Neural network input distributions for isolated track + 2 jet events in the pretag control 
region. The distributions shown are Ai?(lepton-i/m(j2^) (left), PTiW + H) (right). The differences 
in shape are attributable to QCD and are less significant in our higher-purity search regions. 



712 pseudo-data constructed from a sum of background templates. Our expected and observed 

713 limits are shown in Fig. [16] and Table IXVII The limits are expressed as a function of the 

714 Higgs boson mass hypothesis. 

715 The likelihood technique accommodates the uncertainties on our background estimate 
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FIG. 10: Neural network input distributions for isolated track W + 2 jet events in the one SECVTX 
tag region. The distributions shown are Mjj^ (left) and ^^'^^(Loose Jets) (right). The differences 
in the shape are consistent with the uncertainty on our QCD model. 
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FIG. 11: Neural network input distributions for isolated track W+2 jet events in the one SECVTX 
region. The distributions shown are MJ^J^ (left) and Pt Imbalance (right). The differences in the 
shape are consistent with the uncertainty on our QCD model. 



716 by letting the overall background prediction float within Gaussian constraints. We use a 

717 different set of background and signal neural network template shapes for each combination 

718 of lepton type and tag category as a separate channel in the likelihood. We correlate the 

719 systematic uncertainties appropriately across different lepton types and tag categories. 
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FIG. 12: Neural network input distributions for isolated track + 2 jet events in the one SECVTX 
region. The distributions shown are AR(\epton-Umax) (left), Pt{W + H) (right). The differences 
in the shape are consistent with the uncertainty on our QCD model. 
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FIG. 13: Neural Network output distributions for events with one Secvtx tag. The neural network 
output is close to zero for "background- like" events, and close to one for "signal- like" events. The 
open red curve shows the expected distribution of WH Monte Carlo events. The WH expected 
curve is normalized to 50 times the standard model expectation. The plots show isolated track 
events (left) and lepton triggered events (right). 
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Neural Network Output (M=1 1 5) Neural Network Output (M=1 1 5) 

FIG. 14: Neural Network output distributions for events with one secvtx tag and one jet probability 
tag. The neural network output is close to zero for "background-like" events, and close to one for 
"signal-like" events. The open red curve shows the expected distribution of WH Monte Carlo 
events. The WH expected curve is normalized to 50 times the standard model expectation. The 
plots show isolated track events (left) and lepton triggered events (right). 
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FIG. 15: Neural Network output distributions for events with two secvtx tags. The neural network 
output is close to zero for "background- like" events, and close to one for "signal- like" events. The 
open red curve shows the distribution of WH events. The WH curve is normalized to 50 times 
the standard model expectation. The plots show isolated track events (left) and lepton triggered 
events (right). 
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FIG. 16: 95% confidence level upper limit on a{pp — t- WH) ■B{H — t- bb), expressed as a ratio to the 
standard model expectation. The limits were obtained using an integrated luminosity of 2.7 fb~^ 
and analyzing both lepton triggered and I^t + 2 jet triggered events. The dashed line indicates 
the median expected limit. The yellow and green regions encompass the limits in 68% and 95% of 
pseudo-experiments, respectively. The solid line shows the observed limits. 



720 VIII. CONCLUSIONS 



721 Our limit on WH production improves on the previous result by using more integrated 

722 luminosity and extending the lepton identification with isolated tracks. The increase in 

723 luminosity from 1.9 fb^^ to 2.7 fb~^ increases the sensitivity by ~20%. Using isolated track 

724 events provides a ~25% increase in acceptance above the prior analysis. The new isolated 

725 track events combined with minor improvements in background rejection yield a overall 

726 ~15% increase in estimated sensitivity. Our expected limits are expressed as a ratio to the 

727 standard model production rate. The expected limits vary from 4.3 to 43.2 for Higgs masses 
72B from 100 to 150 GeV/c^, respectively. We find no evidence for Higgs production in the data, 
729 and set observed limits at 3.6 to 61.1 for Higgs masses from 100 to 150 GeV/c^, respectively. 
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CDF Run II Preliminary 2.7 fb"^ 
Limits for Combined Lcpton and Tag Categories 
in units of SM cross sections 
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TABLE XVI: Expected and observed limits as a function of Higgs mass for the combined search 
of Tight Lepton and Isotrk events, including all tag categories. The limits are expressed in units 
of Standard Model WH cross sections. 
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